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Abstract 

The energy spectrum and the emission rate of particles emitted from excited nucleus due to 
both the evaporation and the splashing (emission from a cold vibrating nucleus) are calculated. 
We show that the collective motion of the nuclear Fermi liquid is accompanied by direct non- 
statistical emission of nucleons via the dynamical distortion of the Fermi surface. 
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I. INTRODUCTION 



The emission of particles from a collectively excited state of the nucleus as a Fermi liquid 
drop can occur in two ways. First, due to the relaxation processes the collective energy 
is transferred to the intrinsic degrees of freedom with subsequent evaporation of particles. 
On the other hand, a direct non-statistical emission (splashing) of nucleons is also possible 
via the dynamical distortion of the Fermi surface accompanying the collective motion. In 
general the relative contributions of these mechanisms depend upon the magnitude of the 
nuclear friction coefficient. In this work, the limiting cases of the direct (non-statistical) 
particle emission from the non-damped giant multipole resonance (GMR) and the particle 
evaporation from the heated nucleus are compared. The direct particle emission from the 
GMR has been extensively studied within pure quantum mechanical approaches (see, e.g. 

). We exploit here a semiclassical phase-space theory to study the energy 
spectra and the escape width for nucleons directly emitted from a collective state (splash out 
effect). Within this approach, the particle emission is describing in a more qualitative way 
than in microscopic quantum mechanical models. However an advantage of this approach 
is the possibility to take into account the statistical, direct and cascade emissions of the 
particles from the heated nucleus which is simultaneously involved in the collective motion. 
In particular, this approach can be used to describe both the evaporation of neutrons and 
the direct (non-statistical) neutron emission at the nuclear descent from the fission barrier 
to the scission point. Below, we will restrict ourselves to the case of particle emission from 
the giant monopole resonance. A generalization of our approach to the case of an arbitrary 
multipole giant resonances can be done in a straightforward way. 



II. EMISSION OF PARTICLES CAUSED BY THE FERMI-SURFACE DISTOR- 
TION 

For an excited nucleus, the particle emission rate can be calculated from the assumption 
that the nucleons localized in a single-particle mean field of depth Vq and described by the 
phase-space distribution function f{f,'p,t) escape out the nuclear surface if the energy of 
their radial motion exceeds Vq. The fiux density of neutrons emitted from the nuclear surface 
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is given by, see also 

r n rln 

(1) 



surf 



where 6(x) is the function which equals unity for x > and is zero otherwise, gs = 2 is 
the spin degeneracy factor of the emitted neutron, n = f/r is the normal to the surface, 
•fiq = q/q = p/p = flp, V = q/m , m is the particle mass, and E = q^ /2m = p^ /2m — Vq 
is the kinetic energy of the emitted neutrons. We assume that the nucleons are confined in 
a spherical well potential with a depth Vq and time-dependent radius R{t). The integrand 
in Eq. is taken at r = R{t). In the phase-space approach, the distribution function 
f{f,p,t) in Eq. differs from the equilibrium one due to the dynamical distortion of the 
Fermi surface. We use the local-equilibrium approximation and assume the following form 
for the distribution function in nuclear interior jSj 

^' " exp {[{p- mu{f, t)f/2m - /T} + 1' 
where is the chemical potential, T is the temperature and u{f, t) is the velocity field 
associated with the collective excitation of the nucleus. At small temperature T such ex- 
citation corresponds to the zero sound regime in a Fermi liquid 0, Q] within the scaling 



model 



111 ] . The distribution function of Eq. (0) implies the distortion of the Fermi surface 
in momentum space of multipolarities / = and / = 1. This is consistent with the scaling 
model, where the distortions of the Fermi surface of multipolarities / > 2 are absent due to 
the particular choice of the displacement field X {r, t) in the following form 

X {f, t) = a{t) ^0(r), M(f, t) = Qi^ {r, t) (3) 

with 

a{t) = aocos{ujt) (4) 

and ao is the (small) amplitude of oscillations. Considering only monopole density vibra- 
tions, we put 0(r) = r^. The displacement field of Eq. provides the disappearance of the 
highest multipolarities of the Fermi surface distortion with / > 2. Thus, the Fermi surface 
preserves the spherical symmetry under the assumption that the displacement field is given 
byEq. ©. 

Let us relate the time dependence of the chemical potential and the radius R{t) to 
the amplitude a{t). Neglecting the change in the total number of particles A during the 
emission process {A^ 1) we have 



where g = A is the spin-isospin degeneracy factor. Substituting the distribution function 
f{r,p,t) from Eq. ^ into Eq. ^ and using the boundary condition on the velocity radial 
component 



r\r=R{t) 



dRjt) 



we obtain the approximate relations 



(6) 



R{t)^R l + 2a{t) + 2a^{t) , ~ /i 1 - 4 + 8 ^^(t) 



Here, both the equilibrium values R and /i are temperature dependent 
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where Rq and e^? are the corresponding ground state values of the nuclear radius and the 
Fermi energy, respectively. 

Taking into account that the inequality | mu{f,t) Pf is satisfied (see below Fig. 2), 
using Eqs. (j2I) and (jHI) and keeping the lowest non-zero order of the amplitude ao, the basic 
Eq. can be reduced to the following form 



<-^emis(^) 
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7r2 Jo ^ Jo exp {[E + T]{t) - x]/T} + 1 ' 



(9) 



Here, the following notations are used 



Spit) = 2a{t)Ro/vF, r]{t) =Vo~i2{t) ^ X+Aepait), q{E) = 2 ^ep {ep + X + E), (10) 

X = Vo — ep is the neutron separation energy, is the Fermi momentum and Vp = pp/m. 
The energy distribution of the density current of the emitted particles is given from Eq. 
by the following expression 



dJcmisiE, t) 



m 



EQ{E) / dx 



X 



dE n^h^ ' Jo '^'^ exp{[E + 7]{t) -q{E)sp{t)x]/T} + l 

We will also consider the spectral density of emission rate d'^N/dE dt derived as 



'111 



d^N 
dEdt 



AnR 



dE ' 



(12) 
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The derived quantity (fT^ should be averaged over time to smear out rapid oscillations, so 
the observable spectrum has the form 



d?N 



dt 



(13) 



dEdt (27r/a;) 7o dE di 

Here and below the bar means time-averaging. 

In the case of particle emission from the vibrating cold nucleus, the expression (jlip reads 



dJemisjE, t) 

dE 



EQ{E) / dxxQ[q{E)sF{t)x-ri{t)-E]. 



(14) 



We will consider below the particle emission from the giant monopole resonance. In this 
case, the amplitude ao in Eq. (jlj) is small enough to satisfy the condition 77 (t) > (see next 
section). That means that the increase of the Fermi-sphere radius due to the oscillations 
of the nuclear surface is not sufficient to cause the emission of particles, i.e., the particle 
emission appears only because of the dynamic Fermi-surface shift. In this case, Eq. fllljl 
takes the following explicit form 

dJemisi^E y t) 



m 



dE 



E Q{E) Q[spit)] Q[l- e(t, E)] 1 - eit, E) 



where 



V{t) + E 



(15) 



(16) 



Sp{t)q{Ey 

Eq. ()15|) allows us to evaluate the maximal value of the kinetic energy, -Emax, of the 
emitted particles. The result reads 



E ^ 

-'-'max ~ 



-A + 



10 E* 



3 fp 12 A 



\ 5 E* 5 (kpRoy \E* 



(17) 



where kp = pp/fi and E* is the energy of the collective excitation. 



III. NUMERICAL CALCULATIONS AND DISCUSSION 



We will apply Eqs. and (fT3j) to the particle emission from the isoscalar giant 

monopole resonance (ISGMR). The eigenfrequency 00 of the ISGMR in Eq. (jU can be 
found from the classical derivation cu = JC/B, where C and B are the stiffness and the 
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mass coefficient, respectively. Both transport coefficients C and B, can be obtained evaluat- 
ing the collective potential, -Epot, and kinetic, -Ekin, energy of the nucleus within the scaling 
approximation 0. Using Eq. (jS)), we obtain 

^kin = ^rn J d^r p,^u^ ^^AmRl a\t) = ^B a^t) (18) 

and 

B = {12/5)AmRl. (19) 

Here, p^q is the equilibrium particle density. 

The collective potential energy Epot is derived as 

where e is the particle energy density which is related to the static incompressibility K by 

K = 9(p^-^] . (21) 



5p2 



eq 



Using the continuity equation, Sp = — \/ ■ peq X, we obtain from Eqs. Q and 



2 

and 

to = \ 



Epot ~ 2KAa\t) = -Ca^t), C = AKA (22) 



K 



12V (23) 

where (r^) = (3/5)-Rq is the mean square radius of the nucleus. We point out that in the 
case of finite nuclei the incompressibility K is A-dependent. 

Using Eqs. (H)) and (fTHj). one can evaluate the amplitude ao of the Fermi-surface oscil- 
lations. It can be done using the expression for the collective kinetic energy E'kin of the 
monopole vibrations. Averaging Eq. ()18|) over time within the period of oscillations 2tt/uj 
and using the virial theorem we find 

2 - - 1 1 6f E* 

In the case of particle emission from the GMR, we put E* = huo = Egmr, where Egmr is 
the eigenenergy of the ISGMR. In this work we adopt the values of Rq = 1.12 ■ A^/^fm, 
ep = 37 MeV, kp = 1.36 fm'^, Egmr = 82 ■ A-^/'^ MeV and A = 7 MeV. 
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FIG. 1: The spectral density of emission rate fid?N/dE dt for splashing of neutrons from the giant 
monopole resonance (GMR) in cold nucleus (solid curve) and for particle evaporation from thermal 
equilibrated nucleus (dashed curve) versus particle kinetic energy E of emitted neutrons. We have 
used the values of GMR energy Eqmr = 82 • A^'^^'^ MeV, thermal excitation energy Et = Egmr 
and mass number ^ = 90. 



In Fig. 1 the dependence of (PN/dEdt, obtained from Eqs. (fT^ . (fT^ and (fT!?|) . on the 
kinetic energy of the emitted particle E is shown (solid line) for the nucleus with A = 90 
and T = 0. We will compare this result of particle emission from the ISGMR in the cold 
nucleus with the one for the evaporation of neutrons from the thermal equilibrated nucleus 
with the thermal excitation energy Et = Egmr- For the case of thermal evaporation of 
particles, we take fi(t) = X and u{f,t) = in the distribution function f{f,'p,t) of Eq. Q. 
Using Eqs. (|TT|l and (|T^ . we obtain, see also Ref. 



dEdt 



evap 



^ 1 + exp[(i? Iv,- X)/T] ■ ^^^^ 



The dashed curve in Fig. 1 represents the spectral density of evaporation rate (j25|) from 



a heated nucleus with the temperature T = JSEt/A. As seen from Fig. 1, the thermal 
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FIG. 2: The emission rate dN/dt of neutrons from the giant monopole resonance in cold nu- 
cleus with yl = 90 (solid line). The dashed line shows the time dependence of the dimensionless 
parameter, SF{t), of the Fermi surface distortion, see Eq. (|1U() . 

spectrum (dashed line) is much broader than the cold emission spectrum (solid line) and 
has its maximum at a larger energy. 

Fig. 1 shows that the kinetic energy of the emitted particle from a cold nucleus is much 
smaller than the Fermi energy ep, so in ()16|) we can neglect the dependence of q{E) upon 
E. Integrating (fT3jl with q{E) 



Qo = 2 yep {ep + A) over energy E, we obtain an analytical 
expression for the cold particle emission rate dN/dt as 

— = AttR'o Jemis(t) 



^-^^[spit) qore[sp{t)] e[l - Pit)] |l - ^/5(t) + 2p\t) - ^p\t) 



(26) 



where /3(t) = ri{t)/qQSp{t). Fig. 2 shows the emission rate dN/dt of Eq. ()2fi|l for the interval 
of time At < ^-n/io. As seen from Fig. 2 the particle emission from ISGMR occurs as a 
short time splashing. 

The particle emission rate dN/dt can be used to derive the life time, Tcom, with respect 
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to the splashing of neutron from the ISGMR. The quantity ()26|) must be averaged over time 
as in (fT^. resulting in 



For the nucleus with A = 90, one obtains from Eqs. (pUj) and (jTTj) that Xspiash = 6.0 ■ 10~^° 
s. The analogous quantity for particle evaporation, calculated by means of (j25|l . is Tevap = 
5.3 ■ 10^^^ s. The calculations of Xspiash can be improved if the higher multipolarities of the 
Fermi surface distortion, at least the quadrupole, are taken into account, see Appendix. 

IV. CONCLUSIONS 

By use of phase space approach to the particle emission from the nucleus, we have con- 
nected the particle emission rate to the dynamical distortion of the Fermi surface of the cold 
vibrating nucleus. The particle emission occurs here as a classical splashing effect from the 
vibrating liquid drop. We have shown that cold particle emission (splashing) may result in 
the deviation of the observed spectra from the usual statistical ones in the small energy re- 
gion, see Fig. 1. In particular, for the cold particle emission, the emission rate, {dN/dt)sp\s.sh, 
(averaged over time) is significantly large then the corresponding value, {dN/dt)eva.p, for the 
particle evaporation (for the ISGMR ((iA^/(it)spiash exceeds {dN/dt)eva.p by factor about 10). 

The general conclusion of this paper is that the collective motion in a finite nuclear Fermi 
liquid is accompanied by direct (non-statistical) emission of nucleons via the dynamical 
distortion of the Fermi surface. This mechanism of cold particle emission can be also applied 
to the large amplitude motion like a descent of the nucleus from the fission barrier to the 
scission point or to the first stage of the heavy ion collision. 
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APPENDIX A 



We will generalize the local-equilibrium approximation of Eq. Q and assume the follow- 
ing non-equilibrium form for the distribution function in nuclear interior at T = 0, see Ref. 



f{f,p,t) = Q 



(Al) 



where pF{f,t) is the Fermi momentum and v{f,p,t) is associated with the distortion of the 
Fermi-surface in the momentum space. 

To describe the collective excitations with the multipolarity L in the spherical well po- 
tential we will consider h'{r,p,t) as a superposition of plane wave to create a state with a 
good angular moment L, see also Namely, 

.(r.S *) = / da, y„(,4) F(,y . exp|i(* . f- + CO., (A2) 

where 

F{np -Uk) = J2 Yio{np ■ Uk) (A3) 

and z/;(t) = vf''' (e~*'^* + c.c.) is the amplitude of the Fermi-surface distortion of multipolarity 
/ in the momentum space. We point out that the condition / 7^ 1 in Eq. ()A3j) is because 
the shift of the Fermi surface was already extracted in the distribution function f{f,p,t), 
taken in the form of Eq. ()A1|) . We will follow the nuclear fluid dynamic approach, and 
take into account the dynamic Fermi-surface distortion up to multipolarity / = 2 in the 
expansion of Eq. ()A3|) . In this case, the amplitudes z/q and z/2 are related to each other as 
1^2 = a/4/5 z/q, see Ref. jl^. Finally, putting Yl=o,o = l/v^ivr in Eq. ()A2jl . we obtain the 



distortion amplitude i>{f,p,t) for monopole mode in the following form 

u{f,p,t) = uq [jo{kr) - 2]2{kr) P2(x)] e"'^* + c.c, for L = 0, (A4) 

where ji{kr) is the spherical Bessel function, Pi{x) is the Legendre polynomial and x = n-fip. 

The amplitude vq is related to the bulk density variation 5p{f,t). Namely, using Eqs. 
()A1|) and ()A4|) . one obtains the following expression in the nuclear interior 

5p{r,t) = |^-^/(f,p,t)-peq(r) = z/(°)Arojo(A;r)e-^-* + c.c., {r < R, L = 0), (A5) 
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where Nq = 2mpp/'K'^Tt' is the density of single particle states on the Fermi- surface. The 
velocity field u{r, t) can be evaluated from the continuity equation 



Using Eqs. ()A5|) and ()A6|) . one obtains 



— * 



(A6) 



u{r, t) = 2uqNquj 



^^i^fsinM, (L = 0). (A7) 
Po kr 

The wave number k in Eqs. flA4|) . ()A5|) and ()A7|) is derived from the boundary condition. 
The boundary condition can be taken as a condition for the balance, at the nuclear surface, 



between the compressional pressure and the surface tension pressure, see 



15|,|l7|. Taking 



into account the consistent change of the compressional pressure due to the Fermi-surface 



distortion, the boundary condition reads 



id 



Here, z = kR, 



Znjo{Zn) - {fa + ft,)3l{Zn) = 0. 



18 cr 36 /i 



(A8) 



where a is the surface tension coefficient, K' is the dynamic incompressibility given by 



(A9) 



K' = K + K. 



and K is the commonly used static incompressibility 

K = R'' 



,2 



6R^ 



R=Ro 



(AlO) 



(All) 



The additional contribution, K^^, to the incompressibility K' in Eq. ()Ain|l is due to the 
dynamic Fermi-surface distortion effect. The quantity can be evaluated in a gen eral case 
of arbitrary multipolarity of the Fermi-surface distortion and is given by 0, IigI. Il7j| 



12 ep 



fioo(s) 



with n^ofs) 



^1 xPiix)Poix) 

ax 

2 J-i X — s 



(A12) 



where the dimensionless zero-sound velocity s = hu/ppk is found from the Landau's dis- 

nn 

persion relation OQl UJ. In the case of an isotropic Landau's interaction amplitude, i.e.. 



Fo 7^ 0, Fi-^Q = 0, the dispersion relation reads 



^oo(s) = - — . 

-TO 



(A13) 



11 



For realistic nuclear forces, Fq ~ 0, Eqs. ()A12|) and ()A13|) lead to large renormalization of 
the incompressibility with ^ 2K. The giant monopole resonance (GMR) corresponds to 
the lowest solution, > 1; to the secular equation ()A8|) . 

Using Eqs. ()A1|) . ()A4|) and ()A7|) . the flux density Jemis(i) of Eq. ((U) can be evaluated 
beyond the scaling approximation used in Sect. II. 
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